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Abstract 

We study solutions of Ginzburg-Landau-type evolution equations 
(both dissipative and Hamiltonian) with initial data representing col- 
lections of widely-spaced vortices. We show that for long times, the 
solutions continue to describe collections of vortices, and we identify (to 
leading order in the vortex separation) the dynamical system describing 
the motion of the vortex centres {effective dynamics). 



Contents 



Il Introduct 



ion 



2 Ginzburg-Landau preliminaries and results 



2.1 Ginzburg-Landau equation^ 



2.2 Multi- vortex configurations 

2.3 Main resultsi 



3 Proofs! 



3.1 Effective dynamics of vortices: Higgs modell 



3.1.1 Manifold of multi- vortex configuration^ 

3.1.2 Reduced f vortex) Hamiltonian svstemi 



3.1.3 Projections O, 



3.1.4 Splitting! 



3.1.5 Effective dvnamici 

3.1.6 Energy estimates! . 



3.1.7 Approximate conservation of the reduced energy. T-{.(w„) 



3.1.8 A priori momentum bound! 



2 

8 

8 
9 
11 

13 

13 
13 
15 
18 
19 
21 
23 
24 
28 



1 



13.1.9 IterationI 29 

3.2 Effective dynamics of vortices: superconductor mo de] 30 



4.2 Inter- vortex interactio: 



32 

4.1 Approximate static solution propertvl 32 

. . 35 

4.3 Approximate zero-mode property! 38 

4.4 ^ocrcivibL-Ql±he_JIessiai] 39 



4.5 Remainder estimates for GL functional! 41 



^^^endices 



5.1 Appendix 1 



5.2 Appendix 2 



5.3 Appendix 3 



41 

Operators V„\ . . _^ 41 

44 

46 



Proof of Lemma 



Two technical lemma; 



6 References 



48 



1 Introduction 

In this paper we study effective dynamics of magnetic (Abrikosov) vortices in 
a macroscopic model of superconductivity, and of Nielsen-Olesen or Nambu 
strings in the Abelian Higgs model of particle physics. In both cases the 
equilibrium configurations are described by the Ginzburg-Landau equations: 

-A^^ = A(1-J^|2)^ 
curl'^A = Im{ipVA'4') 

where {ip,A) : R"^ ^ C x R'^, Va = - iA, and = V^, the covariant 
derivative and covariant Laplacian, respectively. Equations are the Euler- 
Lagrange equations for the Ginzburg-Landau energy functional 

A) := i jlV^^p + (curlAy + - 1)'^ . (2) 

In the case of superconductivity, the function : ^ C is called the order 
parameter, lipl"^ gives the density of superconducting electrons. The vector field 
A : ^ is the magnetic potential. The r.h.s. of the equation for A is 
the superconducting current. In the case of particle physics, ip and A are the 



2 



Higgs and Abelian gauge (electro-magnetic) fields, respectively. (See jH| for 
reviews, and |JNoj for historical and physics background.) 

In addition to being translationally and rotationally invariant, equations (0) 
are invariant under gauge transformations: 

(^,^)^(e«x^,A + Vx) 

for any % : — > M (solutions are mapped to solutions under this transforma- 
tion). 

We consider various time-dependent versions of the Ginzburg-Landau equa- 
tions (0). The first example is the gradient-flow equations 



= A^v^ + A(i - i^nv' 

dtA = -curt^A + Im{i>\/A^) 



2. , ..^ (3) 



a model in superconductivity theory ([GE[|T]). We will refer to equations (jSj) 
as the superconductor model (they are sometimes called the Gorkov-Eliashberg 
equations or time- dependent Ginzburg-Landau equations). 
The second example is 



dlA = -curt^A + Im{ipVAip) 



2. , (4) 



coupled (covariant) wave equations describing the [/(l)-gauge Higgs model of 
elementary particle physics ( |JTj ) (written here in the temporal gauge). We 
will refer to equations (jH) as the Higgs model (they are sometimes also called 
the Maxwell-Higgs equations). 

The general framework we develop in this paper also applies to coupled 
(complex) Schrodinger and Maxwell equations 

= + A(i - m')^ 

dlA = -curl^A + /m(V^VA^) ^ ^ 

with Rcj > 0, or the Chern-Simons variant of these equations, though the 
implementation for R&f = requires some additional technical steps. 
Finite energy states (■?/', A) are classified by the topological degree 



deg{ip) ■.= deg { — 



\x\=R. 
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where R is sufficiently large (the winding number of ip at infinity). For each 
such state we have the quantization of magnetic fiux: 

B = 27rdeg{ip) E 2nZ, 

where B := curlA is the magnetic field associated with the vector potential A. 

In each case the equations have "radially symmetric" (more precisely equiv- 
ariant) solutions of the form 

^W(a;) = /„(r)e*"^ and ^("^(x) = a„(r)V(ne) , (6) 

where n is an integer and {r,6) are the polar coordinates of a; G M^. As 
r —>■ oo, a„(r) and /n('") converge to 1 exponentially fast with the rates 1 and 
mx := mm{\/2X, 2), respectively: 

fn{r) = 1 + O(e-"^0 and a„(r) = 1 + 0(6"^) . 

At the origin, /„(r) vanishes like r" and a„(r) like r^. Hence 1 — /«(?") and 
1 — a„(r) are well localized near the origin. 

The pair (ip^"'\ A'-"-') is called the n- vortex (magnetic or Abrikosov ( IK \ INoj ) 
in the case of superconductors, and Nielsen- Olesen or Namhu string in the 
particle physics case). Note that deg{ip^'^^) = n. No other static solutions of the 
Ginzburg-Landau equations are rigorously known, though there is a physical 
argument and experimental evidence for the existence of vortex lattices - the 
Abrikosov lattices. 

Observe that (in the present scaling) the length scale for the magnetic 
field (the penetration depth) is 1, and the length scale for the order parameter 
(the coherence length) is where mx = min(v/2A, 2). More precisely, the 
following asymptotics for the field components of the n-vortex were established 
in [P] (see also |.TTj ): as r := |x| oo, 

jW(x) = n(3nKi{\x\)[l + 0(6-™^")] Jx 
S(")(r) = n(3nK,{r)[l - j^, + ©(l/r^)] 

|l-/n(r)| < ce^™-'- ^ ^ 

|/;(r)| <ce-"^^ 

Here j^"-* := /m('?/'(")V^(n)^/'*-"''*) is the ra-vortex supercurrent, and /5„ > is a 
constant. Ki is the modified Bessel function of order 1 of the second kind. 
Since Ki{r) behaves like ce'"^ j ^ for large r, we see that the length scale for 
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j and i? is 1. Note that the two length length scales 1/mx and 1 coincide at 
A = 1/2. Superconductors are referred to as Type / if A < 1/2, and Type //if 
A > 1/2. 

Consider test functions describing several vortices, with the centers at 
points zi, Z2, . . . and with the degrees rii, n2, ... , glued together. An ex- 
ample of such a function can be easily constructed as Vz,x = (V'z.x; ^z,x) with 

m 

^PzA^) = e'^^^^l[ij^^^\x-z,) (8) 
i=i 

and 

m 

Azjx) = ^ A("^)(a; - z,) + Vx(x) , (9) 
i=i 

where z = {zi, Z2, . . .) and x is an arbitrary real function yielding the gauge 
transformation (the integer degrees of the vortices, n = (rii, . . . , rim), are sup- 
pressed in the notation). Define the inter- vortex separation 

R{z) := min \aj — ak\. 

Since vortices are exponentially localized, for large separation /?(z) (compared 
with [min(mA, 1)]~^) such test functions are approximate - but not exact - 
solutions of the stationary Ginzburg-Landau equations. 

When A > 1/2, we take rij = ±1, since vortices with \n\ > 2 are known to 
be unstable f |GSj ). 

Now consider a time-dependent Ginzburg-Landau equation with an initial 
condition fZg.xo ^^^1 ask the following questions: does the solution at time t 
describe well-localized vortices at some locations z = z(t) (and with a gauge 
transformation x = and, if it does, what is the dynamic law of the vortex 
centers z{t) (and of 

We describe here answers to these questions for the superconductor model (jHl) 
and Higgs model (Q. Precise statements (Theorems Q and Ej) are given in Sec- 
tion EH 

Consider the superconductor model © with initial data {ipo: ^o) close to 
some fZo.xo with e~^*^-oV\/^(io) < We show that the solution can be 
written as 

im,m)=vm,xit)+0{e\og'/\l/e)) (10) 
and that the vortex dynamics is governed by the system 

7.,% = -V.W{z) + 0{e'\og'/\l/e)). (11) 
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Here % denotes dzj/dt, iy(z) := ^gl(%x)-E7=i ^^"'^ where := ^gl(V^^"\ ^^"^), 
is the effective energy, and 7„ are the numbers given by 

7n:=^||V^(„)V'(")||^+||cnrM(«)||2. (12) 

In general, these statements hold only as long as the path 'z{t) does not violate 
a condition of large separation: R{z{t)) > log(l/e) + c. In the repulsive case, 
when A > 1/2 and nj = +1 (or rij = —1) for all j, the above statements hold 
for all time t. A precise statement is given in Theorem ^ 

The leading-order term in the r.h.s of ()11|) is of order e (see Lemma fTTl and 
Remark EI). For A > 1/2, the leading order of W{z) for large -R(z) is: 

{const)njnk- 



\/\Zj — Zk\ 



(see Section 

For the Higgs model equations with initial data {ipo, Aq) close to some Wzg.xo 
(and with appropriate initial momenta), we show that 

\\iijit),A{t))-Vzitu(^t)\\H^ + \\{dt^{t),dtA{t)) - dtVzitUit)\\L^ = o(v^) (13) 
with 

ln/z, = -V,W{z{t)) + o{e) (14) 

for times up to (approximately) order log (^) . Here Zj{t) denotes d?Zj{t) / dt^ . 
This result is stated precisely in Section ESI (see Theorem 12)). 

The resulting dynamics of vortices induced by the field dynamics of (z/', A) 
is called the effective dynamics. 

We now outline some previous works on vortex dynamics, including related 
works on the Gross-Pitaevski (or nonlinear Schrodinger) equation 

z^ = -A^ + l(|^p-l)^ (15) 

in a bounded domain, used in the theory of superfluids (see |TTj ) . It is ob- 
tained from (0) by setting 7 = 2 and A = 0. The landmark previous develop- 
ments are summarized in the table below 



6 



Type of 
Eqns 
Type of 
Results 


Superfluid 


Superconductor 


Higgs 


Nonrigorous 


Onsager '49 


Perez- Rubinstein 
'83 (A > 1) 
E '84 (A > i) 


Manton '82 
(A^i) 


Rigorous 


Colliander- 
Jerrard '00 
F.-H. Lin-Xin '00 


Demoulini- 

Stuart '97 
(A^i) 


Stuart '94 
(A^i) 



In more detail, non-rigorous results for the Ginzburg-Landau equation (|TK|) 
without the magnetic component, were obtained by L. Onsager (0), A. Fetter 
(IE]), R. Creswick and M. Morrison (ICM]), J- Neu (|Nj), L.M. Pismen and D. 
Rodriguez (piRj), D. Rodriguez, L.M. Pismen and L. Sirovich (jPRSI), L-M. 
Pismen and J. Rubinstein ( |PiRj ) . N. Ercolani and R. Montgomery ([EM]), 
W. E (jE]), Yu. Ovchinnikov and I.M. Sigal ([OHl). 

Rigorous results are contained in J.E. CoUiander and R.L. Jerrard f |CJj ). 
F.-H. Lin and J. Xin (|L2C|), based on Bethuel, Brezis and Helein (|BBH]). Let 
ip'^ be the solution of Eqn ()15|) with a "low energy" initial condition. Then 
these papers show that as e — 0, the "renormalized" energy density 

_±_(i|vrP + ^(|#P-i)^) 

converges weakly to a sum of (5-functions located at points z(t) := [zi(t), . . . , Zk{t)) 
which solve the Hamiltonian equation i = JVH^z) with appropriate initial 
conditions and Hamiltonian H. Also |CJj prove the Bethuel-Brezis-Helein type 
result Vp > 0, as e — > 

min ||V''-e'"i7zwlki(T2) -^0 

«e[0,27rj 

where is the Bethuel-Brezis-Helein canonical harmonic map with singular- 
ities at zi, . . . , zn and = T^/ U, Bp{zi), and [LXJ show that the rescaled 
linear momentum Im{ijj^'Vipe) converges (on the time-scale 0(1)) to a solution 
of an incompressible Euler equation. The results above describe the dynam- 
ics of the vortex centers, but say nothing about the vortex structure of the 
solutions. 
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In the magnetic case non-rigorous results were obtained in N. Manton ([M]) 
(A PS |), M. Atiyah and N. Hitchin ( \Kn\ ) (A ^ |), L. Perez and J. Rubinstein 
([PRl), and W.E ([Ej). 

Rigorous results were obtained in D. Stuart ([Sj) (A ~ |), and S. Demoulini 
and D. Stuart ffDS]! (X ^ h). 

Final ly, we ment i on the rece nt results |EWl IIWWI (TRFI (TfI IRFI im 

EiEl lEi ESEl EMI EM EMI EB ESI Ell Era IM 

on interface, bubble, spike, and soliton dynamics. 

The rest of the paper is organized as follows. Ginzburg-Landau preliminar- 
ies are given in Sections 12.11 and 12.21 The effective dynamics results described 
above (Theorems Q and |2I) are stated precisely in Section 12.31 Theorem |21 
is proved in Section 13.11 and Theorem ^ in Section 13.21 The key technical 
estimates used in the proofs are themselves proved in Section |3 Technical 
complications are relegated to appendices f Sections 15.11 - 15. 3j ). 

Notation. Here, and in what follows, H"^ denotes the Sobolev space if''(]R^; Cx 
M^) (same for L^, etc.). For /i = (0, a), = (x, P) G L^, (/i, z/) denotes the real 
L^-inner product 



JR2 

Moreover, we will use the same symbol to denote the real inner-product in 
X L^: for ^ = (^i,6), V = ('71,^2), we write 



L^-norms are denoted with a subscript p: || • ||p = || ■ The letter c will 
denote a generic constant, independent of any small parameters present, which 
may change from line to line. 

2 Ginzburg-Landau preliminaries and results 
2.1 Ginzburg-Landau equations 

The Ginzburg-Landau energy functional Sql (see Q) is a smooth functional 
on the following affine space of configurations of degree n: 




(16) 



(17) 



:= {(V-, A) : ^ C X | (^^ a) - (V-^"), A^")) G H^} 
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where {ip^'^\ A^"''^) is the exact n- vortex solution of the Ginzburg-Landau equa- 
tions (see (j2I))- The variational derivative £'qj^{iI),A) is the (negative of the) 
right hand side of the Ginzburg-Landau evolution equations Q (or Q). 

With the notation u = {ip, A), the superconductor model equations can 
be written as 

dMt) = -ScMt)). 

We consider solutions of © satisfying u = [tp, A) e C^(]R+; X^")) (see jHE] for 
existence theory). 

It is convenient to write the Higgs model equations (@]) as a first-order 
Hamiltonian system. Introduce the momenta 

{n{t),E{t)) := {-dMt),-dtA{t)) 

{E{t) is the electric field). The Hamiltonian is 

nitl^,A,7,,E):=£GL{^,A) + l [ {\7t\' + \E\'}, (18) 

a smooth functional on the space X*^") x L^. The space X^"^ x L^, viewed as 
a real space, admits the non-degenerate symplectic form 

a;(e,r/) = (e,r^) (19) 

where (■, ■) is the real inner product on the tangent space to X defined 
in (fTTjl . and J is the symplectic operator 




(in block notation). Setting w := (■?/', A, vr, i?), the Higgs model (jH) is equivalent 
to the equation 

dtw{t) =m'iwit)). (20) 

We consider solutions in the space w E C^(M;X^"^ x L^) which conserve the 
Hamiltonian functional Ti. (see |BMj for existence theory). 

2.2 Multi-vortex configurations 

We begin by constructing a manifold of multi-vortex configurations, made up 
of collections of widely-spaced vortices "glued" together. Such a collection is 
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determined by m G Z+ vortex locations, z = {zi, . . . , z^) G IR^*" and m vortex 
degrees, n = [rii, . . . , rim) G Z™, associated with these locations (the latter will 
often be suppressed in the notation), together with a gauge transformation. So 
the manifold we construct may be parameterized by a subset of M^™ x {gauge 
transformations}. 

Recall ('(p^^\ A^"'^) denotes the equivariant, ra-vortex static solution of the 
Ginzburg-Landau equations (see ©)• To a triple z G M^™, n G Z™, x : ^ 
M, we associate the function 



where 



and 



Vz,x-= iipz,x,Az,x), (21) 



i=i 

Here (ni, . . . , rim) are the fixed topological degrees of the vortices, rij G Z\{0}. 
For given z G R^™, the gauge transformations will be of the form 

m 

i=i 

with X £ i?^(M^; M). The gauge transformation is taken to be of this form to 
ensure that vz,x lies in X^"'\ 

Given a vortex configuration z = (zi, . . . , Zm), the inter- vortex distance is 
defined to be 

R{z) '■= min \zj — Zk\- 

l<j<k<m 

To ensure that our multi- vortex configurations are approximate solutions of the 
Ginzburg-Landau equations, the inter-vortex separation will be taken large. 

In the Higgs model case, momenta must be included. To do this, we first 
introduce the "almost zero-modes" . Define the gauge "almost zero-modes" 

Gf'^^:={^,d^)vz,x (22) 
for 7 : ^ M, and the gauge-invariant translational "almost zero-modes" 

T;?'^^ := id... + - ^.), d^))vz,x. (23) 
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From explicit expressions for Gf'''^ and Tjf^^ (see dH and (03), one can 
deduce that G^^'^\Tj^'^^ G H"^, provided 7 G H"^^^. Then for momentum 
parameters 2 = {Pi, ■ ■ ■ ,Pm) ^ M^™ and C ^ i/^(lR^; M), we define the (vr, E) 
(momentum) component to be 

771 

^z,,,u,--=Y.PrTt-"'' + Gf'^'eL\ (24) 
i=i 

We will often denote the full set of parameters by cr := (z, x,R, C) and 0z,x,p,c 
by (f)a. 

An important role will be played by the interaction energy of a multi-vortex 
configuration (see Section IT^ : 

m 

where, recall, E^") := ^Jl^ £gl(^^"\ ^^''^)- Due to the gauge invariance of 
Sgl, this interaction energy is independent of the gauge transformation x- 

2.3 Main results 

The main result in the superconductor model case is as follows: 

Theorem 1 Suppose A > 1/2 and rij = +1 (or rij = —1) for j = 1, . . . ,m. 
There are do, d,eo > such that for < e < eo the following holds: let 
{ilj{t), A{t)) solve (0) with initial data satisfying 

\\{ijo,Ao) -vz^^^Jhi < doelog^/\l/e) 
with e-^(^oVV^Uo) < doe. Then for t > 0, 

\\im,m)-vmxit)\\H^ < de\og'/\l/e) 
for a path vz(t),x(t) ^ Mas satisfying 

buM^) + V.,Vr(^)| < de' log3/4(l/e), (26) 

m 

Wdtxit) -Y,z,{t) ■ A("^)(- - z.mim < de'\og'/\l/e). 
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Here 7„ is a positive constant, given explicitly in (jl2j) . 

For the Higgs model equations, we have the following result: 

Theorem 2 Suppose A > 1/2 and rij = +1 (or rij = —1) for j = 1, . . . , m. 

Let a(e) be a function satisfying < a{e) « log"^^^(l/e). There are 
do,d,T,eo > such that for < e < eo, the following holds: let w(t) = 
{ip{t), A{t),Tr(t), E(t)) solve ^EW> '"^^^^ initial data satisfying 

\\ii^o,A) -Vz^^xoWm + ll(7ro,^o) -0io,xo,£o.foll2 < doe\og^^\l/e) 

with e-«(^o)//R(^ + |£J2 + < doe. Then forO<t<^^ log (^) , 

\\{m,m)-Vmait)\\m + \Mt)^m)-M)\\2 < da{e)V~eW/\l/e) (27) 
for a path a{t) = CW) satisfying, for all j , 

\h - Pil + \ln,Vo + V,^.iyU)| < e«(e) W/'(l/e) = o(e) 

- E ■ ^^"''^(^ - ^^■) - ^11^^^ + ll^*Clk-= < log^/'(l/e) = o(e) 

(28) 

/or any s > 0. 

Remark 1 The inter-vortex force is of size e: 'WW{z) = 0(e) (see Lemma[T7\ 
and Remark{^ . 

Remark 2 The condition A > 1/2 and rij = +1 in Theorems^ and\^ en- 
sures that the inter-vortex interaction is repulsive, and therefore that the inter- 
vortex separation does not become too small in the given time interval. In 
fact the theorems apply, without these restrictions, for any initial vortex con- 
figuration whose evolution (namely i28\} or \2b]) ) preserves an appropriate 
large-separation condition. In the Type II case (X > 1/2), this condition is 
e-^(^i^)) / ^ R{^z{t) < e (and \2{t) \ + ||C(t)||i2 < t in the Higgs model case). In 
the Type-I case (X < 1/2 this condition must be appropriately modified (see 
Remark\^ of Section \J~^ ). 

Remark 3 In Theorem\^ since \zj\ < c^e over the time interval < t < 
T = log (^-^^ ; vortices can move a distance 

y^T = rlog^ ~ R{z{0)) » 1. 
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3 Proofs 

We start by proving Theorem |21 for the Higgs model in Section ITTl The proof 
is considerably more involved than that of Theorem ^ for the superconductor 
model. The latter proof is sketched in Section 13 .21 

The proof of Theorem |21 given in the following section is based on a se- 
ries of propositions and lemmas. Propositions summarize our geometric 
construction, and Lemmas [T]|7l whose proofs are left to Section HI provide the 
(elementary) analytic building blocks. Several technical lemmas are relegated 
to appendices. 

3.1 Effective dynamics of vortices: Higgs model 

In this section we prove Theorem|21 Let w{t) solve (j^UI) with w E C^(K; X*^") x 
L^). In what follows, we denote 

X := H\R^; C x M^) x ^^(R^; C x R^). 

3.1.1 Manifold of multi- vortex configurations 

We begin by defining the manifold of multi-vortex configurations. Let 

S := {(z,X,E,C) I z G M2",x-z-^z e i/2(M2;M),pG M2"^,C e H\R^;R)}, 

where z-Az '■= SjLi ' Aj with Aj(x) := A^'^^\x — Zj). This set is a manifold 
under the explicit parametrization map 5 : 1^2,1 S defined by 

'5: (z,X,2,C)^^(z,X + z-Az,l),C)- (29) 

Here 

Yr^s ■■= ffi^™ X //^(R2; R) X M^"^ X H'{R'^; R). 
We define an open domain in S by 

:= {(z,x,2,C) e S I e-^®//R(i) < e, |p| + \\C\\m < V~e}. 

For each a := (z, x, 12, C) ^ ^! introduce the multi-vortex configuration 

:=(z;z,x,0<x)eX(")xL2 (30) 
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(recall V'z^^ and ^o- are defined in 1)211) and ()24|1 ). Finally, we define the space 
Mmv := {w, k G C X 

The map 7 : ^ X*^"^ x given by 7 : cr ^ i^o-, parameterizes Mmt,, so that 
^rav = li^e)- It is easy to check that 7 is C^. It is shown in Section I^TTl that 
for all cr e Se, its Frechet derivative D'~f{a) : Tg-Se ^ X is one-to-one. Hence 
Mmv is a manifold. 

For each a G S^, the tangent space to Mmv at Wo- will be denoted by 
Tw^Mmv It can be identified with a subspace of X; specifically, T^^Mmv = 

For use in computations and estimates below, we introduce convenient 
bases in T^T,^ and T^^M^v In terms of the coordinates in ()29p . the basis in 
To-Sg is given by 

{dz^, + {zi ■ d^^Ai, d^),d^(^), dp^^,d^(^)}. (31) 

We denote the coordinates of a' G T^T,^ in this basis by cr'coord ^ ^2,1- Define 
the map : ^2,1 — ^ by 

^^coord ■■= D^{a)a'. (32) 
For cr(t) a path in S^, this definition implies 

where cT(t) is the coordinate representation of the vector dtcr(t) G To-(j)Se: 

d(t) :=(z(t),(fWxW,EW,5*CW), 
with z(t) := dz{t)/dt, p(t) := d2.{t)/dt, and 

m 

i=i 

Let (9^, := (9^-^ + The basis for T^^M^v (which is the image of 

the basis (|^T|) under D^{a)) is given by: 

:= df^^'Wa, rf- := dp^^w^, := d^^^)W^, := (33) 

14 



Note that the tangent vector t^ - is defined by differentiating Wa- covariantly. 
The point here is that {d^)"^w„ lies in x for any m, while dzW^^ does 
not. Exphcit expressions for these tangent vectors are given in ()105|) - ()108|) . 

For a vector a G M^™, we will set a ■ r* := a^rf^ for r* = r^^-, r'f-, and 
for a function 7, set (7, r*) := J 'y{x)T*dx for r* = r^, r|. As a result of 
these definitions, and the relation 

dt = z-dz + {dtx,dj^) +:g_-d^+ {dtC,d^) 

= i-di+< cf d^>+^-d^+ {d,C, 5c), ^^^^ 

we have 

^^coord = z' ■ + (x', r>^) + p' ■ + (C, r«), (35) 

where ^^^^^ = (z',x',e',C') e >2,i. 

In what follows, all of our computations are done in these bases, and we 
omit the subscript "coord" from the coordinate representation cr'^oord ^ vector 



3.1.2 Reduced (vortex) Hamiltonian system 

As was discussed above, the Maxwell-Higgs equations constitute a Hamilto- 
nian system on the phase-space X*^") x with Hamiltonian psp . Our goal 
below is to project this Hamiltonian system onto the manifold M^v (more 
precisely, onto TM^v) with the smallest error possible. Below we describe an 
equivalent Hamiltonian structure on the parameter space 1^2,1 which is used in 
our analysis. We begin by setting 

Xr,s := H'^iM.'^; C X M2) X H'iR'^; C x M^) 

(note that X = Xi^). The operator To- has adjoint A^- (with respect to the 
l^2m X X M^™ X inner-product on Yr^s, and the real x inner-product 
on Xr^s) given by 

A^:^^ {D^w^,0 (36) 
or, in our coordinates in Tg-Se, 

A. : e ^ i{r:,,0, {r^,0, {r!,,0, (r^O)- (37) 

It is shown in Section 15.11 that To- and A^- are bounded uniformly in a G 
between the following spaces: 

To- : ^r,s A^r-l,s-l; (38) 
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for any r and s satisfying min{r, 1) > s — 1, and 

^o- '■ ^r,s ^r-l,s-l; (39) 

for any r and s satisfying min{s, 1) > r — 1. (In the Physics hterature the 
operators To- and Ag- are called bra and ket vectors, with notation To- = \Dwcr) 
and Act = {Dwa\-) 

Define the operators 

K := AJ-'T^ : Yr,s ^ ^-2,^-2, (40) 

where s — 1 < min{r, 2). 

Relation (pn|) shows that V* = —Va in the sense of the 1? inner product. 
The operators define a symplectic form on 1^2,1 by 

u,,la' ,a"){G) :=((7',Ka"). 

The non-degeneracy of this symplectic form follows from: 

Proposition 1 (non-degeneracy of reduction) For e sufficiently small, and 
a G Se, the operator is invertible. 

Proof: The invertibility of the operator for sufficiently small e follows from 
the following expression, shown in Section 15.11 

t ) <«' 

where 

r._( D 0{e\og'/\l/e)) \ ( 0{^e) \ 

""-yOieW/^l/e)) K j' ""'-y-Oi^ey J' 



/ 0(6logV2(l/e)) \ 

''^-Uo(elogV2(l/e))* J' 

0(ye) stands for an operator whose norm is bounded by c-y/e, D is a matrix 
of the form Dj^^im = Inj^ji^km + 0(elog^^^(l/e)), and K is the operator K := 
—A + |V'z,xP- Since D and K are invertible, the operators Ua and Va are 
obviously invertible if e is sufficiently small. □ 
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The symplectic form uo^edW ) cr") and the reduced (vortex) Hamiltonian 
h{a) := T-C{Wfy) give a reduced Hamiltonian system on 1^2,1- The corresponding 
Hamiltonian equation is 

& = V-^Dh{a). 

This equation will turn out to be the leading-order equation for the dynamics 
of the parameters a. The next proposition computes the Hamiltonian h{a) 
exphcitly. 

Proposition 2 If e^^^-^ ^J R{z) < e, then 

m ^ m ^ 

j=i j=i (42) 

+ 0(elogi/^(l/e)(bP+||C||^)) 

where, recall, := £gl(^^"\ ^^"^); In = ||| V^(n)V'^''^||2 + \\curlA'-''^\\l and 
W{z) is defined in / t^,5|) (it is computed to leading order in Section \J~^ . Fur- 
ther, we have 

Dh{a) := D,H{w„) = {VW{z),Q,r2,K()+0{eW'\l/e){\p\ + \\C\\2)), (43) 
where 7 ■ p denotes (71P1, . . . , 7mPm) ■ 

Proof. We begin with auxiliary computations establishing the approximate 
orthogonality of the tangent vectors introduced above. To this end, we record 
explicit expressions for T^-f'^^ and G^^'^\ which follow readily from defini- 
tions (EH) and (1221): 

^ _(e^x[-Q^(nO(^ _ ^,)](V^,^)('^^)(x - z,), B^-^\x - z,)ei) (44) 

and 

G(2-x) = (,^^^^^^ V7). (45) 

Here i?*-"^ = V x A^""^ is the n-vortex magnetic field, := (0, 1) and 65 : = 
(-1,0). 

By the above explicit expressions, the exponential decay estimates (|7j), and 
Lemma IT^ we see 

KT}.^'^\Tif>^^)l<C6l0gV2(l/6) 
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when j ^ k. When j = k, we compute 

+ {BjJer,BjJes) +0(e), 

and the leading term is easily computed to be 6rs where 7„^. is given in (fT^ . 
Thus we have the approximate orthogonality relation 

D,,,im := (T,|'"\t£'"^) = InM^^rn + O (e log^/^ ( 1/e) ) . (46) 
A similar computation yields 

|(i;?'"\Gf>^))|<celog^/^(l/e)||7||2. (47) 

Finally, the corresponding relation for the approximate gauge modes (see (j221)) 
is 

(Gf Gf^) = (7, (-A + li^zJ'X), (48) 

a straightforward calculation. 

Now using Wo- = {vz,x, 'Pa) (with (p^ defined in (j211)), and 

110.11^ = P,.P..(T,r\ T,|-^) + {Gf'^\ Gf-) - 2p,,{Tji'-\ Gf-)) 
together with we obtain ^ and □ 

3.1.3 Projections 

Here we construct operators used to engineer a convenient splitting of (j2Uj) . 
We define the operator : X ^ T^^Mmv as 

Qo := ToV-'Aor\ (49) 

Due to the expression for Va- in (HOI), we see that is a projection, = Qo-, 
and it satisfies 

KerQo = (JT^,M„,)^ (50) 

and 

q: = -jQ.jj. (51) 

Finally, we list two estimates which follow readily from the definitions 
above: 

||g.||x.,.^x.,. <c (52) 
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for any r and s satisfying s < min{r + 1, 1), and for a = cr{t) a path in Sq-, 

\\[Qa,dt]\\x^x < c\\&\\yi_o (53) 

where, recall, a = (z, df-x, R, dtC) ■ To obtain this estimate one uses rela- 
tion dni. 

3.1.4 Splitting 

The next proposition establishes a coordinate system (adapted to the projec- 
tion Qo-) on a tubular neighbourhood of Mmv Let, for < do < 1, 

S° := {(z,X,E,C) e S I e-«®//R(i) < rfoe, N + \\C\\m < doV~e} 
(which parameterizes a manifold somewhat smaller than Mmv)- Set 
Us := {w e X^") X I \\w - w„\\x < S, for some a G 

Proposition 3 (coordinates) For e sufficiently small, there is 6 » e, and 
a map 

satisfying Qs{w)(yU! — ws{w)) = for w & Us. Moreover, DS{w) is bounded 
uniformly in w eUs- 

Proof: The proof is an application of the implicit function theorem. Define 

g-.Us^iTfl-^ r_i,o 

by 

g{w, a) := A^J"^(u; - w^). 

One can check that this is a map. Obviously, giw^, cr) = 0. Note that, due 
to dnil), D„g{w^, a) : ^2,1 ^ Y-ifi is given by 

D^g{wa,cr) = -A^J^^r^ = -14. 

which is invertible for cr G with e sufficiently small. So the implicit function 
theorem applies to provide a map w i— S{w) from an if^-ball of size 5 of a 
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given Wo- G M^v into S, satisfying g{w, S{w)) = 0. Allowing a to vary in 
we can construct such a ball about any such w^- 

Using the definitions (jH^ and (jHBj) of the operators and Ao-, the explicit 
expressions (jlUSj) - (jllOj) for the basis {t"" ,t^,tP,t''}, and expression (HH), one 
can check the following: there is Sq independent of a G such that for all 
w G BxiWfj] (5o), the norms 

llK"^ll>'-l,0^^2,l' I|r<7||y2,l--^1,0' ||Aa||xo,l^>'-l,0' 

||-D<jA^||xo,ixy2,i^y-i,o' \\DcT^AY2,ixY2,i^Xifi, ll-D^Aa||y2,ixy2,ixXo,i-^y_i,o 

are bounded uniformly in a. This fact implies that the balls on which the 
maps 5* are defined can be taken to be of uniform size 5 « -^/e, which implies 
S{w) G Eg. Thus we obtain a well-defined map S : w y-^ S{w) from the 
tubular neighborhood Us into S^, with g{w,ws{w)) = 0. This map obviously 
satisfies also = 0. The uniform boundedness oiDS{w) follows 

readily from the formula DS{w) = —[Da-g{w, S{w))]~^Dwg{w, S{w)) and the 
uniform estimates mentioned above. □ 

Now suppose w(t) solves the Higgs model equations with initial data 
w{0) = wq as specified in Theorem El In particular, we have w{0) G Us- Let 
< Ti < oo be the time of first exit of w{t) from Us- For < t < Ti we may 
write 

wit) = w^^t) + m (54) 

with Wa(t) G Mmv, and Qa{t)C,(t) = (by choosing a(t) = S{w{t))). By our 
choice of initial data, 

\\m\\x<cUo\\x<cdoe\og'/\e), (55) 

where '■= w{0) — Wa-o- Indeed, using and the equation w{0) = w^^ + 
we find 

^(0) = - w^(o) + ^0- (56) 

Next, since a{0) = S{w{0)) and ctq = S^Wa^) (see Proposition EI) , and since 
w{0) — = C,o, Proposition El gives 

Iko-o-(0)||y2,, <c||^o||x. 
The last estimate, together with the estimate HZ^o-Wo-Hy^ -^^x < c implies that 

which, together with yields (jSHI)- 
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3.1.5 Effective dynamics 

Insert the decomposition into the equations (p])) and expand in a Taylor 
series to obtain 

dtw^ + dt^ = + Lai + N„{.i)] (57) 

where L^. := 7i"{wa) is the Hessian of H at w^, and 

N^iO := n\wa + - ^'K) - Lai 

consists of the terms nonhnear in ^. Apply the projection Q„ to (jFfj) and use 

Qadtw^ = dtw„ (since dtw„ E T^^Mmv) to obtain 

dtw^ - Qjn'iw^) = QalSL^i - dti + IN^iO]- (58) 

This equation governs the effective dynamics of the parameters o"(t). The 
terms of leading order are on the left hand side. We now show, starting with 
the nonlinear term, that the right hand side is of lower order. 

Lemma 1 (nonlinear estimate 1) For a E Se, and ^ := (^1,^2) £ xL^, 

NAO - ( ('^"'J*^ ) 

with 

\\iNMii)\\H-^<csmirH^ + Ui\\u 

for any s > 0. 

This lemma is proved in Section f4.5[ From now on, we fix s > (s = 1/2, 
say). Thus by we have 

\\QjN,mL^><H~^<cim\i+urx)- 

To minimize writing in the rest of this section, we make the additional 
assumption 

mix < 1, (59) 
which we shall justify later. Then the above estimate becomes 

\\QJNamL^xH-s<cU\\j,. (60) 

Using the fact that Qfj^ = and the bound we have 

\\Qadti\\x< c\\a\\y, JiWx. (61) 

To bound the remaining term on the right hand side of ()58|) . we need the 
following lemma whose proof is given in Section 14.31 
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Lemma 2 (approximate zero-modes) For cr G and any /5 G x L^, 

we have 

\\LaQal3\\L^^L^ < Cv^||/3|U2xL2. (62) 

Fix 1] ^ L"^ X L'^. Using tlie symmetry of L^, and ()(j2j] . we have 
\{V,QJLM\ = \{L^QJv,0\ < U\\mxL4LaQJv\\H-^xL^ 

< Cy/e\\^\\xM\L'^xL2 

and hence 

WQJLMIl^xl^ < cVeUllx. (63) 

Collecting (jHOl), and (jHS)), we obtain a bound on the right hand side of 
the effective dynamics law (j3Hj) : 

\\dtw^ - Qjn'{w^)\\L2,H-s < c{V~e + UWx + ||cr||n,o)||e||x. (64) 

Finally, we translate (jMjl into parametric form, in order to remove a from 
the r.h.s., and to see that it yields in the leading order. For a G C^(R; E^), 
we recall 

dtw„ = (65) 
where a = {z,p,d^x,dtC)- Next, using (j49|l . we find 

Qjn'iw^) = T^V~'AM{w^). (66) 

Now the definition of A^., (jHSl), implies that 

A^n'iw^) = D,n{w^) = dM^)- (67) 

The last two equations yield 

Qjn'iw^) = T^V-'DMwa). (68) 

Comparing with (jUHj) . we obtain 

dtw^ - Qjn\w^) = T,{& - r-'5Xa)). 

Now the mapping properties and (jHIIjl . and the fact that Ag-Fo- : yi,i-s — > 
is invertible, imply that 

||cr - V7^5^/i(o-)||yi,i_, < c\\dtW^ - QJH'{w^)\\l2^h-''- 
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Since Dh{a) = {H! {wa), D^Wa) = 0(e^/^) (this estimate is part of Lemma 
below), (jnH) implies that 

||a - V-'DM^)\\y...-s < c{V~e+ Uhmh. 

Now using and ()4H1 . we arrive at 

V-'DM^) = feC,-7-' ■%W^,0) + O(elog^/2(l/e)(|n| + \\Ch)), 

where we have used the notation 7"-*^ • dzW = {•j'^dz^W, . . . , 'Jn^dz^^W). Com- 
bine this relation with the estimate above to obtain finally 

m m 

\m-vAt)\ + E ip.w +7n7v.,H^(z(t))i + Wxit) - at)\\m + mmm 

<c[{V~e+m\x)m\x + e^"W\l/e)]. 

(69) 

3.1.6 Energy estimates 

Our remaining task is to control the remainder ^{t) for long times. The idea 
is similar to techniques used to prove orbital stability of solitary waves in 
Hamiltonian systems (see, eg, jWl IGSSj ): exploit conservation of energy - in 
this case both for the PDE ()20|) and for the leading order effective dynamics (PH|) 
- in order to control the fluctuations. We begin with a Taylor expansion of 
the Hamiltonian: 

n{w^ + = n^w^) + («;.), + l^O + ^^(0 (7o) 

(this equation deflnes i?o-(0)- The following lemma, proved in Section \4A\ 
allows us to control ^ by the Hamiltonian. 

Lemma 3 (coercivity) For e sufficiently small, cr G S^, and ^ G kerQ^, 

-m\x<{i,L^^)<cm\\. 

c 

Using this lemma, together with conservation of the Hamiltonian, in (jTOJ, 
we obtain 

m\l < c[n{w{o)) - niw^) - {n'iw^),o - R.m 
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The following lemma bounds the super-quadratic terms on the right-hand side 

of dZD). 

Lemma 4 (nonlinear estimate 2) For a G S, 

\RAo\<c{urx+urx)- 

This is proved is Section 1^31 To control the terms linear in ^, we need another 
key lemma: 

Lemma 5 (approximate solution properties) For a G S^, we have 

1. \\H'{wa)\\HixL2 < cv^e 

2. m'iw^Mn^ = \\S'advz,^)\\m<celog'/\l/e) 

3. \\Qjn'iw^)\\L2xL^ < celog'/\l/e), where := 1 - Q,. 

This lemma is proved in Section I^^Tl Using the third statement of the lemma 
and (jHU, we find 

l(7^>.),e)i = i(H'K),g.OI 

= |(JQ.J^'K),0I 

<\W'H{w^)\W^LA\i\\x ^ ^ 

<ce\oil\\le)\\iU. 

Collecting Lemma El Lemma El and (j72p in (jTlj) (and remembering the 
intermediate assumption < 1), we obtain 

llell^ < c[7^(«;.(o)) - n{w^) + (6 W/^(l/e) + \\if^m\x 
+ (elog^/^l/e) + 11^(0) lU) 11^(0) lU]. 

3.1.7 Approximate conservation of the reduced energy, T-C{w„) 

It remains to control li.{wa{Q))—l-i{wa)- The estimate below involves a delicate 
estimate of the contribution of the nonlinear terms. 

Proposition 4 Let M(t) := suPq<^<( II^('5)|U- Then 

\niv^io))-niv^^t))\ < ctv^M(t)(eW/2(l/e) +M(t)) + ViM^(t). (74) 
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Proof: First, we differentiate in time and use the effective dynamics law (|58|): 



d 
dt 



(75) 



where we have used the fact that (Qo-J)* = —QaS- 

We start by estimating the first inner-product on the right-hand side. First 
we exploit the fact - Lemma El part 2 - that the first component of T-C'{wcr) is 
smaller than the second: using QadtC, = [Qa,dt]^, and (jHSI), we have 



To deal with the second component, we have to exploit a key cancellation. 
This is expressed in the following lemma, which can be considered a refinement 
of both dnH) and of 

Lemma 6 For a G C"'^(M; S^) and Q^-^ = 0, we have 



This lemma is proved in Section l5^ Combining the above estimates with (j69|) 
yields 



|(H'K),Q.[JL.e - m)\ < c[e'/'W/\l/e) + eUWxmWx. (76) 



Finally, we must control the second inner-product on the right hand side 
of (j75|) . This is problematic since, so far, we have control over (the second 
component of) the first factor only in (see Lemma Ej), and over the second 
factor only in H~'^ for s > (see Lemma Q). The solution is to isolate the 
worst term and to use the detailed structure of the equations to deal with it. 

First we claim that T-C'{wcr) is of the form 



\{[n'iw^)]u[QA^L^^-dtO]i)\ < cehg 



(lA)llellx(v^+||a||y,J. 




+ mh+^re{\z-2\ + \\d^x-c\\Hm\\x. 



rest 



(77) 
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with := {(,d^w^), and H'^^g^ satisfying the estimate 

WKeJH^-'xH^ < cV~e. (78) 
Indeed, fl77j) - (fTSjl is easily obtained from the exphcit expression 

n'iw,) = i£'aL{v^,^)A.), (79) 

where (p^ = PjkTjf'^^ + Gf"^\, and estimates \\Sgl{'"?l,x)\\h^ ^ ce log^/'^(l/e) 
(Lemma part 2) and \p\ + ^ v^- Thus 

{n'{w„),QjN^{0) = -{Q4-H^ + mLtlN^iO) 

= {Gf''\ {NMii)) - iQJKest.NAO)- 
So by (UHl), (EOI), and (jHI), we have 

\{n'iw.),QJN^iO) - {Gf ''\ {NMm < cVml- (80) 
Next, we single out the worst term in the nonlinearity: 

where 

II -^rest II xL^ < client- 

Recall here that we are writing ^ = (6,6), 6 = HiiiiA)-, and ^2 = (^7r,^£:)- 
Hence 

KGoiv.e.*)i<c|icik.|ieiil. 

Then in light of (j8(J|) . we have 

m\w„),QJN„{i)) - {VC,Im{^^VA^J^))\ < cVeUWx- (81) 

It remains to estimate (V ■ (, Im{^^V (note that the estimates 
available to us so far control the first factor in (and no better) and just 
fail to control the second in L^). The key is to recognize this quantity as 
(essentially) a time derivative. Using the basic equation (jHTjl and (fTHjl compute 

j^{CJ'm{l^i^)) = {dtC, Irn{^^^^)) + (C, /m(^^[-940^]^ + [^gl(^z,x)]i^ 
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We estimate each term on the RHS as follows: 

\{CJm{^{N^),]^))\<c\\C\\Hm\'x + mx) 
\{C,Im{^4S'advz,dW)\ < celog'/\l/emnAmx 
\{C,Im{^^d,[<P^l))\ < c(|e| + + ||9*C||2)(|z| + \\cfxh) 

+ \\dtCh)\\C\\mm\x 

and 

|(C,/m(e.[9,^z,x-[0a]^]))l 

= KC,Me.[fe-l)..[T,|'^^]^ + [gI^^ >]))! 

<9rx-C 

<c(iE-z| + ii^x-ciii/Oiicik^iieiu. 

We write 
with 

\{CJrniC^^Lt))\<c\\C\\mmx- 
Collecting these estimates yields 

< cm\x{m\\m-Amx + VelUWx + 6log^/^(l/e) + IeI 

+ V~e{\i\ + WS^xh) + \R-i\ + \\S^X- CWm]) 
<cv^llellx(||e||x + 6log^/^l/e)), 

using ()69p . Noting that 

and combining (|7H|) . (fTnj) . (jHH), and (jH21), we have 

||[7^K) + (C,/m(e"^e.))]| < cV-emUUWx + e\og'/\l/e)). 
Integrating this in time, and defining M(t) := suPq<^<( leads to 
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Returning to (f7H|) . we obtain 



WmWx < c[tV~eM{t){M{t)+e\og'/\l/e)) + M{t){e\og'/\l/e) + M'{t)) 
+ v^M2(t)M(0)(elog^/^(l/e)+M(0))]. 



It now follows that there is a constant r' > 0, such that for < t < 



In particular, the intermediate assumption ()59|1 is justified. 
3.1.8 A priori momentum bound 

We wish to iterate the above argument to extend the time interval. The 
problem is that the vortex velocities can, in principle, grow to size et » -^e 
if t >> l/\/e (this would mean we leave the manifold Mmv, and many of the 
above estimates fail). We show here that this does not happen. To this end 
we use the approximate conservation of the reduced (vortex) energy H{wa), 
together with the repulsivity of the interaction energy. Indeed, since we are in 
the "repulsive" case (A > 1/2 and rii = ■ ■ ■ = rim = ±1), we have the following 
lemma, which is proved in Section 

Lemma 7 (interaction energy) For R{z) large, 



Here Cjk > are constants. 

Remark 4 One can see from this expression that like-signed vortices repel, 
while opposite- signed vortices attract. 

Conservation of energy for the PDE, T-L{w{t)) = TC{w{0)), together with 
the decomposition (jK^ and a Taylor expansion yields 



mm 





(84) 



+omrx+\\mrx)- 
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We saw above (in ^) that \{n'{w^),0\ < ceW/^(l/e)||^|U, which gives 
n{w^) - 7^K(o)) < cmWx + ||e(0)||x)(elog^/^l/6) + UWx + ll^(O)llx)). 
By estimates (j42j) and K > c for some c > 0, we have, for e sufficiently small, 

m 

for some constant a > 0. In light of Lemma [7J the assumptions on the initial 
conditions in Theorem El imply ^.{10^(0)) — Ylf=i E^"'^^ < c'd^e, for some c'. We 
choose do < where a' is a constant to be chosen below. So provided 

ciiimh + ||e(0)||x)(elog^/^l/e) + + \\m\\x)) < ^6, (85) 

we have ap-C^Wa-) — Yl'jLi E^'^^^] < a'e, and therefore 

W{z) + \i,\' + \\C\\l^<a'e. (86) 

So by (jHll), if a' < min(cjA;), then as long as condition holds, |pp + 
IICIII-i < e, and i?(z)e"-^® < e. Hence a e S^, and e Mmv 

In particular, this estimate shows that Ti > r'/ ^/e. Hence, we have shown: 



Lemma 8 There are r' > and d > 0, such that inequality ^EB) holds for 
<t < t' / y/e, provided 

\\m\\x+\\m\\x<dV~e. (87) 

3.1.9 Iteration 

We may iterate Lemma |H1 for as long as the conditions a G S^, and ||^(t)||x < 
d^/e hold. Iterating N times starting with ^(0) and satisfying \\C,{0) \\x < doy/e 
yields 

WmWx < Cc^elogi/2(l/e) for < t < r'N/^e, 

where C is another constant. The condition (jH7|) limiting the number of iter- 
ations, ensures both that (jH3|) holds (so that a G remains true), and that 
the remainder in the effective dynamics law is sub-leading order. Thus we can 



29 



take = a(e)/ye for any aQ^/e < a{e) « log (with oq > 1). This 

gives a total time interval of length 

where r = r'/ logc, over which we have the bound 

\m\\x<Ca{e)V~e\og'/'il/e) (88) 

for < t < T. 

Finally, equation (jHnjl implies 

m-m + \m+T' ■ Vzw^(z)i + ncfxw - mwm-^ + mmH^ 

<cV~emx<cea{e)log'/'{l/e) = o{e). 
This completes the proof of Theorem |21 □ 

3.2 Effective dynamics of vortices: superconductor model 

Here we just sketch the proof of Theorem [T] since it proceeds as above. The 
important difference is that we can control the remainder for all times. 

The set-up is as follows. For the superconductor model our manifold of 
multi-vortex configurations is taken to be 

= {vz,^ I e-^®//R(i) < e, X e ff|(M'; M)}. 
A solution u{t) = (tplt), A{t)) of (jS)) is decomposed as 

U{t) = Vz{t),x{t) + ^(t) 

with -Pz,x^ = O5 where Pz,x denotes the orthogonal projection from onto 
the tangent space T^^^Mmv Substituting this into ^ yields 

where Lz,^ := ^g'l('^z,x)- "^^^ equation governing the effective dynamics of 
z(i) and x(t) is derived by applying the projection P^ ,^ to this: 
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To estimate the RHS, we have the following properties 

||Pz,x^z,xellL^<celogi/2(l/e)||^||Hi 

\\Pz,^N^zJO\\H-s<cm\'m + Urm) 

\\Pz,Aam<ci\i\ + \\^x\\Lm\\m- 
Combining these yields 

Hvz,^ + Pz,>cS'gl{vz,>c)\\h-s < c(6 log^/2(l/e) + |z| + \\cfxh + mmMWrn 

which in parametric form reads (recall the notation ji := {'jniZi, . . . , 'Jnm^m)) 

|7z + V.l^(z)| + WdhWm-' < c{e\og'/'{l/e) + mmMWH^. (89) 

In order to control for all time we need 

Lemma 9 There is 6 > such that for R{z) sufficiently large, 

This lemma is proved in Section 14.41 

We use the fact that the main part of the energy difference, i^(fz,x + ~ 
^(%,x)' namely L^^xC) ^ is a decaying quantity. Compute 

+om + \\dah)m?H~)- 

Now we use 

K^'(%x),%xe)l<ceW/^(l/e)||e|ki, 

i(iv„^jo,^z,xe)i<ciieii^.(iieik^ + iieii^o, 

and 

together with Lemma El to obtain 

(| + ^)(e,%,xe)<iieiii.[c(iieiki + 11^11^1 

+ |z| + ||5a||2)-5/2] + celog^/^(l/e)||^||H.. 
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(90) 



we have 




Finally using 



c{^,L^,^0<mm<-{^,L^,. 



0, 



we find 



M\t) < c[e-(^/'^)*M2(0) + elog^/'(l/e)M(t)] 
and so if M(0) = 0(e log^/^(l/e)) we have 

\m\\m<celog'/\l/e) 

for all t, as long as and e"-^®/ V-^fe) < e hold. By dHHl), we see 

|z + Vzl^(z)| + WdhWm^s < ce'\og'/\l/e) 

so the intermediate assumption is justified. Finally, in the repulsive case, 
£{u) - ^^"L^ < ce implies e'^^^VV^fe) < e holds for all t. □ 

4 Key properties 

In this section we prove the lemmas used in the proofs of Theorems ^ and |21 

4.1 Approximate static solution property 

Proof of Lemma O The main fact we use here is that since we consider 
the Type-II regime (A > 1/2), the effects of the magnetic field and current 
dominate those of the order parameter at large distances. 

In what follows, a subindex k will denote an equivariant field component, of 
degree rifc, centred at Zk. eg, := V^^'-'H" -^fc), (V^V')fc = - 
Zfc), etc. 

We first prove 
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Lemma 10 



\\£'GL{vz,^)h<ce-^^^^/R"\z) 



(91) 



Proof: The proof is a computation using the fact that m'^"^^ = A'^'^^^) 
satisfies the Ginzburg-Landau equations, together with the exponential de- 
cay ((7j). We start with 

[^'gl{vz,x% = -^Az,^^z,x + A(|^z,xr - l)^z,x- 
Using gauge covariance and the covariant product rule, we find 



3«X 



^Az,x^z,x = e 
A little computation plus ((Tj) yields 



. (92) 



(n/i)-i-E(/i'-i)is'E^- 

j=X j=l j^k 



mx{\x-Zj\+\x-Zk\) 



(93) 



Using (|^ and (jHSJ, together with the fact that m*^"^^ solves the Ginzburg- 
Landau equations, we arrive at 

I[^gl(^z,x)]^(x) - [E^^'%,{x)\ < c^e-^(l^-^^l+l^-^'=l) 



where 

^(z,x) _^.x ^( -Q ^,)(V^^), ■ (V^^), 

Using Lemma [T21 (in Appendix 3, Section 15. 3^ with a = P = 2m x > 2, 
7 = 5 = 0, we obtain 



< cg-m,R(z)^^^^3/2 e-^®//R(i)- (94) 



We turn now to 



[^gl(^z,x)]a = curlBz,x - jz, 



Observing that curlBz,^ = "^^^iCurlBj, and 
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using the Ginzburg-Landau equation curlBj — jj = 0, invoking an equation 
similar to (jHSl) for Ylk^j fl using ((Zj), we arrive at 

\[S'gl{Vz,x)]a{x) - Ef''\x)\ < C ^-n.^(\.-z,\+\x-z,\)-\x-z,\ 

j,k,i distinct 

where 

Estimating as above gives 

II[^gl(^z,x)]a - Ef^\\2 < ce-^^©i?(z)3/2 e^^^/v^- (95) 
Using (|7j) again, we obtain the foUowing pointwise estimate for E^-'^^ = 



|E(^-)|<cJ:- 



(1 + |X 1)1/2 (l + |x-Z,|)l/2- 

Applying Lemma IT21 with a = (3 = 2 and 7 = 5 = 1 yields 

||^(z,x)||2 < ce-^^^^ /R{zf/\ (96) 

Then ((111)-® yield (El). □ 

Now we consider the manifold of approximate solutions for the Higgs model 
equations. Parts 1 and 2 of Lemma El follow immediately from the expression 

(cf. m) 

together with Lemma ITUl and the fact that a G implies |p| + \\(\\h^ < e and 
e-^^y^/R{zj < e, which implies e--^®/^(z)^^^ < celog^/^(l/e). The refined 
statement, part 3 of Lemma follows from the fact that for a G S^, 



where Pz,x denotes the orthogonal projection onto the span ofT^f^^ and G\ 
(see Eqn. (ITTK|) of Appendix 2) and so, since Pz,x[PjkT^f^'^ + Cf'^^] = 0, 

\\Qjn'iw,)\\L2^L^<ce\og'/\l/e) 
as required. This completes the proof of Lemma El □ 



34 



4.2 Inter-vortex interaction 

The reduced energy is a function of the vortex positions alone: 

m 

In this section, we compute - to leading order in the vortex separation - Vr(z), 
and Viy(z), the inter- vortex force entering the effective vortex dynamic laws. 

Proof of Lemma [7| Noting that B = curlA, we re-write the Ginzburg- 
Landau energy as 

^GL{A,^P) = l [ {\VA^P\' + B' + ^m'-m. 

For (■?/', A) = {4'z,x^ ^z.x)) we have 

m 

j=i k^j 

and B = Bi + ■ ■ ■ + Bm- So plugging into £gl and using the notation ji : = 
Im{'4}iV Aiipi) and fi := we find 

m 

^gl(^z,x, %,x) = ^^"'^ + + 



where 



and 



^ 771 



i=l fc^i JT^/ 

+ 1 / E(/i - - 1) + E (/i - - - 1) + ■ ■ ■ ]• 
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For each term in Rem, the integrand is bounded by e~*^™™(™^'^)l^~^jl+'"^l^~^'^l) or 
Q-{mx\x-Zk\+\^-zj\+\x-zi\) ^ g^^fj after integration, is << e~"^®/A/-R(z) (using 
Lemma IT^ and mx > !)• Using the Ginzburg-Landau equation curlB = j, we 
can re-write the leading-order term as 



A computation gives (-A + l)B = n{2{l - a)ff + a'(l - p))/r > 0. By ©, 
|(-A + l)B{x)\ < ce-"^^!^-!, and = CnUe-^^^ / + 0(l/|x|)], c„ > 0. 

Applying Lemma El yields 

LO = i ^Qn.n.^iL^ / e-(^-^'=)/l^-^'=l(2(l - a,) Mi, + <(! - A'))/rrfa; 

\/\Zl~Zk\JR2 



+ 



e 



-m) 



Lemma [7| follows. □. 

Now we turn to the estimate of the force: 



Lemma 11 We have 



V.,Wiz) = Y.^,n,C,:^^==^^^ + o(e-«(^)/v^) (97) 

7^ ^/\Zj - Zl\\Zj - Zi\ 

as R{z) —>■ oo. Here Cji > are constants. 

Proof: By the definition of W{z), V,,„W^(z) = (^^^(t^z.x), t/I'^'V Equa- 
tions (jni) and (jnni) imply that 

V,„l^(z) = T.g'^)) + o(e-^-®/ v^), (98) 

where is defined in the proof of Lemma ITUl We first compute 

n(Z,x) r^(Z,x)i \ _ ik 



where (recall the notation ipkix) := il)'^"'^\x — Zk), etc.) 



jk 
O^lm 
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First, we note that aj^ = af;^. Second, we use ((Tj) to conclude that if / 7^ j and 
I 7^ k, then {af^l << e~-^®/ a/ -R(z)- It remains to compute a^-^. We rewrite 
to get 



■jk 
(X ■ 



and use 



s r^j,k 



JJ^ j2 _ ^ ^g-mA-max(|x-2j|,|x-Zfe|) 



to conclude that aj^ = aj^ + o(e -^©Z ^jR{z)), where 



Writing everything out in terms of the vortex profiles fj and Qj and taking the 
real part, we find (applying Lemma IT^ again) that 



/ 








r 


k - 



n{l-a)ff' 



+o(e-^®//R(i))- 



Now using the fact that [xrn{Jx)s — {Jx)mXs] equals if s = m, — 1 if (s, m) 
(1,2), and 1 if (s,m) = (2, 1), and summing over s, we arrive at 



~ik 



n{l — a) 



J k 



n{l-a)ff 



[X 



T,)„ + o(e-^®//R(i))- 



Now we apply (a slight variant of) Lemma IT^ to obtain 

Q-\Zj-Zk\ 



\/\Zo-Zk\ 

+ o(e-^®//R(i))- 



(^T^fc)™ / e^-^'^-'^^l\-'^-'\\-a)frirdx 



K2 



Thus 



k^l 



-\zi-Zk\ 



\Zl - Zk\ 
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:(ir^fc)^ + o(e-^®//R(i))- (99) 



The computation of (-Sf''^, [^/m '^^]^) similar, but simpler. We just report the 
result: 



^-\ZI-Zk\ 



{E^\ [Tt^'^U) = 5^4n,^^=(^T^,)^ + o(e-^®/v^)- (100) 



Combining ^ and (fTUnD with (jHED yields □ 

Remark 5 Similar computations can be made for the Type-I case, A < 1/2. 
In this case, 

W{z) = 0{e-"''^^^'^) 
as R{z) oo, and the inter-vortex forces are attractive. 

4.3 Approximate zero-mode property 

Proof of Lemma El Set Lz,x '■= ^Gik^lL.x) ■ -^^^ ^^-^ write 

-^z,x = ^5 + ^0) 

where Lj := ^gl(£/xo)«^"'H- - ^i)). X(j) ■= X + 'Ek^j^i' - ^k), and Vq) is a 
multiplication operator satisfying 

\V(j){x)\ < ce-™"^-^^l^-^'=l. 

The notation g^u stands for the result of acting on m by a gauge transformation 
7. Recall the translational modes T^'^^ are given in (jl^)- Using the fact that 

we get the easy estimates || LjTj^''^'' || 2 < ce~^®, and 

\\y{j)^jk h See ' '. 

Thus 

ll^z,xT,|'"^l|2<c6W/^(l/e). (101) 

To deal with the gauge modes, G^'^^ := {l^d^vz,^, we use f)13|) . which 
gives 

%,,Gf^) = (^7[^^GL(^Z,x)]V',0) 
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and so 

\\L^,^Gf^%<ce\oi'\l/e)hh. (102) 

Now by (jS21), and piop . RanQa = T^^Mmv consists of vectors of the 

form (a ■ T^'^^ + Gf'^\ 0^2(7^)) with a G M^^ ^nd e H\ This, together 
with 

and (fTTITIl and (fTU^ . yields Lemma El □ 
4.4 Coercivity of the Hessian 

Proof of Lemma El Suppose i] := {jj^, rjAj is orthogonal to each approximate 
translational zero-mode, Tj^'^\ and to the approximate gauge zero-modes, 
G^'^^ (which means Im{ilJz,xVip) = V ■ ?7a by an integration by parts). Set 
L := Lz^x- fi^^^ S°^^ show 

{v,Lv) > cillr^ll^i. 

Let {Xj} be a partition of unity associated to the vortex centres. That is, 
SjLo ^1 ~ supported in a ball of fixed radius about Zj (j = 1, . . . , m), 

and xo is supported away from all the vortices. By the IMS formula ( |('FKSj ). 

L = J2xjLx,-2Y,\'^Xj\'- 

We can choose {xj} such that |Vxj| < cR~^, where R := -R(z). As in Sec- 
tion jOl set 

and write, for each l<j<m, L = Lj + V(j). Since 

l%(^)l<c$^e-l— '=1, 

we can choose {xj} so that || V(j)Xj||oo < Cy/e, and so 

{xjV,Lxjv) > ixjV^LjXjv) -cV^WvWl- 

for 1 < j < m. Also, since xo is supported away from all the vortices, 

{XoV,Lxov) > C2\\xov\\h^ 
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for some C2 > 0. Thus 

m 

{v,Lr]) > '^{XjV,LjXjV)+C2\\xoV\\m-c{V^+R'^)\\v\\m- 

Now let {Tjk} {k = 1, 2) be the exact translational zero-eigenfunctions of Lj 
(see |GSj for a discussion). We have 



\{Tjk,XjV) \ < ce, 

and 

So by the n- vortex stabihty result of IC^Sj (for rij = ±1), we have 

{XjV,LjXjV) > csWXjVWm - ce\\T]\\l, 

and so 

{r],Lr]) > [c^-ciVe + R-^)]M\Hi > cMlm (103) 

for e sufficiently small. 

For the Higgs model, the linearized operator acts as the identity on the 
momentum components: 



L.:=7^>.)=(Y I)- 



Observing that Qo-^ = implies -Pz.x^i = 0{^/e) (see Eqn. pi5|l in Appendix 
2), we have 

{^,L^O>b-o{V~mMm + \Ml 

This proves Lemma |21 (the upper bound is straightforward). □ 

Proof of Lemma ini Set L := i^z.x- First observe, using LPz,x = o(l), that 

Now since + AW^i-^l^ < c, for any < 5 < 1 we have 

m,Lo = 5m,Lo + ii~5)m,Lo 

> 6{A^, AO - csmi. + (1 - 5)(ci - o(i))iieii^. > m\i^ 

for 6 and e sufficiently small. □ 
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4.5 Remainder estimates for GL functional 
Proof of Lemma |^ For v = {ijj, A), set 

Here 

and Sqi^{v) is the Hessian of Sql at v (which we don't write out exphcitly 
here). After some computation, we find, for ^ = {r],a), 



+ Ai?e(^r7)|r/|2 + ^I^Hr^p + ^jr^l^j 

and so using Holder's inequahty, and the Sobolev embedding \\g\\p < Cp\\g\\Hi 
in two dimensions, we obtain easily 

\Rvm<c{mm + mm). m) 

□ 

Proof of Lemma 12 The most problematic term in Ny{C,) is of the form .^V^, 
so we will just bound this one (the rest are straightforward): 

II^V^IIj:,-. = sup |(r/,eVOI<sup||r/e||2||Ve||2 
ll»?ll«''=i 

<csnp\\vmumm<curm 

where 1/p + 1/q = 1/2 and q is taken large enough so that C L^. □ 



5 Appendices 

5.1 Appendix 1: Operators Va- 
in this section we consider the key operators := A„^^^T„, where T„ and A^. 
are given in ()35|1 and ()37|1 . and we show for them the relation PT|1 . implying, in 
particular, the invertibility of V^. We also prove the auxiliary properties (p?H|) 
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and (jH^ of the operators T„ and A„. To this end, we use the following explicit 
expressions for the basis vectors for the tangent space T^^M„ 



'■mv • 



r!, = ^^''\SJ^, (105) 



^j^ = (0,T,ra (106) 

r^ = {G'i'''\Fl), (107) 

r^. = iO,Gff^) (108) 

where 

SJk ■■= m,, + tA,,)[<l,^Ud.J<l>^]E) (109) 

and 

Fl :=a^(,)</),= (z4[0,]^,O). (110) 

In what follows we omit the super- and sub-indices (z, x) and a. Using 
Equations p05|) - piO|) it is not difficult to verify properties (jHHj) and (jH^ . For 
example, to show (jH^ we calculate using definitions ^7\i and ()105p - (jllOp . 

where ^ = (^1,^2) 

7(x) :=(G,,,ei) + (F5,,6), 

Clearly |a| < c||.^||xrs '"5-^5 ^"^^ similarly for j3. Furthermore, due 

to (03), 

and due to (jllO|) . 

Recall that := E7=iPi " ^j"'''^ + '""^ and that C e H\ Using that ^/^(R^) . 
H'{M.'^) C i/'^'(M2) for r' < min{s, 1), we obtain 

hll/z'-i < c(||^vIIh'-i + II^aIIh'- + II^ttIIhO < c||^||h'-xh», 
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provided r — 1 < min{s, 1). Similarly, we have 

Ihll//"-! < c||6IIh- < c||^||h«"!/xh»- 
Summing this up, we conclude that 

provided r — 1 < min{s, 1), which implies ()39p. ()38p is obtained similarly. 

Now we use explicit expressions ()105|) - ()108|) for the basis ()33p in order to 
establish equation (jH}. Equation (PT|) follows from the relations (r, J~^r) = 0, 
where r = r^^-, r?^-, r^, or r|, and the relations 

{r^.J-'r^i) = -{^J-'rQ = InAkh + 0(e log^/^(l/e)), (111) 
(^^•,r\^), (r^,rV,^,) = O(v^), (112) 

= (ri,rV,^,) = 0(elogV4(l/6)), 
-V) = -(rf,J-V>^) = -;r., (114) 



r; 



|2 



where i^'^y is the integral kernel of the operator K = —A + \4'z,x\ 

We will not present here proofs of all the relations (jlll|) - (jll4|) . but rather 
illustrate our arguments by establishing two of the relations, say {Tfj,3~^T^) 
and (Tj^-, JI~^Ty) (see flll2|) and pi3|) ). In what follows, we omit the super- 
scripts in T^'^\ Gf^'^\ Sji^, and Fj^^, and the subscripts in t/^z.x and (pa- Using 



Equations ()l()5p and p()7|l . we obtain 

{r^',J-'r^) = {T,k,FsJ-{S,,,GsJ. 

Using the explicit expressions (jUj), (jlSI), (!lU9j) . and (jllUj) for the vectors on 
the r.h.s, we compute 
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where il:^jk){x) = [Yli^j 'ip^'^^^x — ~ ^j)- Recalling the defini- 

tion of (j), and using \p\ + \\C\\h'^ < we conclude that pi2p is true. 
To prove that {T^jjI^^r^) = 0(e), use Equations (|1U5|) and (|l(J7p to obtain 

Now using (jUj) and and the equation curW'^'^^ = /m('?/;'^"^V^(n)'i/'*^"^) and 
estimate ((Tj), we find 

- curlB^''^\x - Zj)\ 
= |/m[(n - zi)\^ - mVA\k^Y^^Kx-z,)ij^^^\x - zj)]\ 

and using Lemma fT2| we see \{Tjk,Gsy)\ < celog^^'^(l/e). 
5.2 Appendix 2: Proof of Lemma [Bl 

We first note that using Qo-^ = and therefore Qadti = [Qaidt\i, and us- 
ing and (UHl), we obtain 

(7^'K),g.9,o = (0.,([Q.,5*]O2) + o(eW/^(i/6)||eilxll^lln,o) 

Using the above expressions for V^, To-, and Ao-, and differentiating Qo- with 
respect to t leads to 

(-A+|V.z,J2)-i((V[9.,,+ ^ ^^'^ ^ \] + {drX,dx))Gf'''\i2) 

+om\\xm + mh + Mi\ + \\^x\\m))). 

Using 0, = j9,tTif^^ + gJ-'^^ we find 

{MQ^^dtih) = (v5;-i^^ ,6)+0(v^||eiU(|sl + l|5*CI|2+v^(|z| + ||ffx||^^i))) 
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and so conclude 

|(7i'K),Q.a,0-(v^;-i^;z ,6)1 < c(v^||e|U(|E| + ||5*C||2+v^log^/'(l/e)(|z| + ||9^xbO))- 

We turn now to computation of {H'^Wa-), Qa^L„^). To do this, we have to 
refine the above computations, and compute Qo- up to 0(e). We find 

where Pz,x denotes the orthogonal projection onto the span of the vectors 
T^f''^-' and G^'^\ and Qu and (521 are ©(i/e). We will need the exphcit form 
of Q21: 

Q21 = D-\T^^'^\ ■)S'^ + (-A + \^PzJ')-'{G^-'''\ ■)^" 
- D-\S-, ■)T(^'^) - (-A + lij^J'r'iF'^, ■)G,. 

We have 

= ([g.(0.,-^'(t;z,J)]i,£:"(t;z,x)ei) + ([g.0.,-^'(^z,x))]2,6) 
= (Pz,x0 + 0(6=^/^ log^/'(l/e)), S"{v^,^)i,) 

+ (-Pz,x^'(t;z,,) + g2i0 + 0{e'/' log^/2(l/e)), 6) 
= (£"(t;z,x)^z,x0,6) - (^'(^z,x),%xe2) + (Q210,6) 

+ 0(6^/2 logV2(i/e)||e|U). 

Now use ^"(fz,x)%x = C'(elog^/^(l/e)) and the fact that 

= Q.^=(Pz,xl + O(v^))e 
which implies Pz,x^j = 0(v^||^||x) to find 

{n'iw^),QJL^O = (g2i0,6) + 0(e3/Mogi/2(l/e)||e|U). 
From the form of Q21 given above, and the fact that {T'-^'^\^2) and (G^^'^^^) 
are 0{^/e\\C,\\x), we see 

(H'K),QaJi^.O = (/^-'(T(^'^\0.)5'^+ (-A + |z^z,xr)"'(G(^'^\0<x)i^",6) 
+ 0(e^/2W/2(l/6)||e||x) 

Combining this with the above computation of {TC'{wa), QadtCl, with the facts 
that IeI = |p+V,l^(z)|+0(e) and \z\ + \\^x\\m = \i-Mdtx-C\\m+0{^t) 
proves Lemma El □ 
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5.3 Appendix 3: Two technical lemmas 

Lemma 12 Let < a < p and < 5,'j < 3/2. Then 

^a\A^-p\x-a\ e-"l"l / |a|-i/2 a = l3 

Proof: We prove only the case a = (3, since the remaining cases follow from 
Lemma IT^ below. Define / by 

r p-a\x\-f3\x-a\ 

J^2 \x\^\x-a\^ 

oo T i'2tt -ar-pJr^ + \a\^-2r\a\ cos(6») 



/o Jo (^^ + - 2r|a| cos(e))'5/2 ' 

Changing variables to m = 1 — cos(6') and t = r/\a\, and using a = /?, we 
estimate 

dt du e-^m+^(.t-i)-+2tu) 



+ / =: Ji + /2 + /s 



/ < 



'0 Jo J1/2J0 Ji Jo 
Using {t - If + 2tu>[l-t + tu]"^ (for < m < 1), we have 

io ^'^ Jo 



(117) 



ce-H/^r^re-^Vv/^c.. ^^^^^ 



t^-'/' Jo 



a . 



Now estimating ^/ {t — 1)^ + 2tu > a/(1 — t)^ + m for t > 1/2 and changing 
the variables of integration as x = |a|(l — t) and y = \a\^/u we obtain 

Jo Jo + y)' 

rW\ r\a\ -ay^/(2^x^+y^) 

< ce-^l^'l/lap-^ / dx / dy- 
Jo Jo 



(3,2 + ^2)5/2 



2\a\ r2. 
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We have 

^5lg-"»-cos2(e)/2 ^ ^ / "■^ p-r(l+s)a/4 



Hence 



Finally, 



h < 



h < c^=. (119) 



(r + 1)7-1 Jo ^ (r2 + 2n)^/2 
2^..|.| , dr ri ^^,-.\a\ir^V^) (120) 



/o (r + 1)7-1 (^2 + 2t;2)5/2 

< ce-°l"l/|a|'"'^. 



Estimates (fTTTIl - imply (fTTI)l) . □ 



Lemma 13 Suppose b{x) is a function satisfying \b{x)\ < ce~™'^' f 
m > 1, and e{x) is a bounded function with asymptotic behaviour e{x 
Cie"!""!/ a/|x|(1 + 0(l/|x|)) as \x\ oo. Fix z eR"^ and set 



or some 



I{z) := / b{x)e{x — z)dx. 

Then 

I{z) = cie-l"l/v^ / e"-"/l"l6(x)rfx[l + 0(1/1^1)] (121) 

as \z\ oo. 

Proof: Choose a with 1/m < a < 1. Let Da\z\ denote the disk of radius a\z\ 
about the origin. We have 

I / 6(x)e(x - z)dx\ <c e-'^Wdr < c|z|e-'"'"l^l « e'l^l/kr (122) 

iR2\D^|,| Ja\z\ 
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for any p. On Da\z\i \x\ « \z\^ and we Taylor expand: 

\z-xt^l^ = \z\-^'^\l + 0{\x\/\z\)l 
e-\'--\ = e-\'\e''-'/\'\[l + 0{\x\^/\z\)\, 

yields 

I h{x)e{x - z)dx = e~\'\ / ^/^ I e'-''\'\h{x)[l + 0{\x\^ /\z\)]. (123) 
Estimates (fT^ and (fnH|l together yield (fTT[|l . □ 
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